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ABSTRACT 



Aims. We present the results of the study of the substructure and galaxy content of ten rich clusters of galaxies in three different 
superclusters of the Sloan Great Wall, the richest nearby system of galaxies (hereafter SGW). 

Methods. We determine the substructure in clusters using the 'Mclust' package from the 'R' statistical environment and analyse 
their galaxy content with information about colours and morphological types of galaxies. We analyse the distribution of the peculiar 
velocities of galaxies in clusters and calculate the peculiar velocity of the first ranked galaxy. 

Results. We show that five clusters in our sample have more than one component; in some clusters the different components also 
have different galaxy content. In other clusters there are distinct components in the distribution of the peculiar velocities of galaxies. 
We find that in some clusters with substructure the peculiar velocities of the first ranked galaxies are high. All clusters in our sample 
host luminous red galaxies; in eight clusters their number exceeds ten. Luminous red galaxies can be found both in the central areas 
of clusters and in the outskirts, some of them have high peculiar velocities. About 1/3 of the red galaxies in clusters are spirals. The 
scatter of colours of red ellipticals is in most clusters larger than that of red spirals. The fraction of red galaxies in rich clusters in the 
cores of the richest superclusters is larger than the fraction of red galaxies in other very rich clusters in the SGW. 
Conclusions. The presence of substructure in rich clusters, signs of possible mergers and infall, and the high peculiar velocities of 
the first ranked galaxies suggest that the clusters in our sample are not yet virialized. We present merger trees of dark matter haloes in 
an N-body simulation to demonstrate the formation of present-day dark matter haloes via multiple mergers during their evolution. In 
simulated dark matter haloes we find a substructure similar to that in observed clusters. 

Key words. Cosmology: large-scale structure of the Universe; Galaxies: clusters: general 



1. Introduction 

Most galaxies in the Universe are located in groups and clusters 
of galaxies, which in turn reside in larger systems - in superclus- 
ters of galaxies or in filaments crossing underdense regions be- 
tween superclusters. According to the Cold Dark Matter model, 
groups and clusters of galaxies form hierarchically through the 
merging of smaller systems (Loeb 2008; Knebe & Miiller 2000, 
and references therein). The timescale of the evolutio n of groups 
depends on their global environment (Tempel et al. 2009). As a 
result, the properties of groups depend on the environment where 
they are embedded, richer and more luminous groups are located 
in a higher-den sity environmen t than poor, le ss luminous groups 
dEinasto et alj l2003allbl l2Q05h iBerlind et all |2006|) . An under- 
standing of the properties and evolutionary state of groups and 
clusters of galaxies in different environments, and of the proper- 
ties of galaxies in them is important for the study of groups and 
clusters, as well as for the study of the properties and evolution 
of galaxies and larger structures - superclusters of galaxies. 

Detailed knowledge of properties of clusters of galaxies is 
also needed for comparison of observations with N-body calcu- 
lations of the formation and evolution of cosmic structures. 

Galaxies and galaxy systems form because of initial density 
perturbations of different scales. Perturbations of a scale of about 



100 h~ l Mpdj give rise to the largest superclusters. The largest 
and richest superclusters, which may contain several tens of rich 
(Abell) clusters, are the largest coherent systems in the Universe 
with characteristic dimensions of up to 100 h Mpc. At large 
scales a dynamical evolution takes place at a slower rate, and the 
richest superclusters have retained the memory of the initial con- 
ditions of their formation, and of the early evolution of structure 
dKofman et al.ll 987). Rich supercluste rs have high density cores 
that are absent in poor superclusters dEinasto et alJl20 07b). The 
core regions of the richest superclusters may contain merging X- 
ray clusters (Rose et al. 2002; Bardelli et al. 2000; Belsol e et al.1 
2004). The formation of rich superclusters had to begin earlier 
than that of the smaller struct ures; they are the sites of early 
star and galaxy formation (e.g. iMobasher et al] 2005b. and the 
first places where systems of galaxies form (e.g. lVenemans et all 
2004; Ouc hi et alj|2005h . 

Among the richest galaxy systems, the Sloan Great Wall, 
which is the richest s ystem of galaxie s in the nearby Uni verse 
(IVoeelev et all 120041; iGott et all 120051 iNichol et all 12006). de- 
serves special attention. The SGW consists of several rich and 
poor superclusters connected by lower density filaments of 
galaxies, representing a variety of global environments from the 
high density core of the richest supercluster in the SGW to a 
lower density poor supercluster at the edge of the SGW. The su- 
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perclusters in the SGW differ in morphology and galaxy content, 
which suggests that their formation and evolution has been dif- 
ferent (Einasto et al. 2010, in preparation, hereafter E10). 

Our aim in the present paper is to study the properties of 
the richest clusters of galaxies in the SGW. Our clusters come 
from all the superclusters in the SGW so that we can compare 
the dynamical state and galaxy content of the richest clusters in 
different environments, from the high-density core of the richest 
supercluster to clusters in a poor supercluster in the extension of 
the SGW. 

The present-day dynamical state of clusters of galaxies de- 
pends on their formation history. One indicator of the former or 
ongoing mergers betwee n groups and clusters is the presence of 
substructure in clusters dKnebe & Mlillerl 12000. and references 
therein). The substructure affects the estimates of several clus- 
ter characteristics, the dynamical mass and mass-to-light ratio 
among others (see Biviano et al. 200q, for a review). Mergers 
also affect the properties of galaxies in clusters. The well-known 
morphology density relation tells that early type, red galaxies 
are located in clusters (in the central areas) while late type, blue 
galaxies can preferentially b e found outside of rich clusters, or 
in the outskirts of clusters flE inasto et al. 1974; Dre sslerlll980t 
Butc her & Oemlerll978HEm"a sto & Einasto fl987h . An old ques- 
tion is whether the properties of galaxies depend on the clus- 
tercentric radius or on the local density of galaxies in clusters, 
or on both ( Whitmore & Gilmore 1991; Huertas-Companv et al. 
120091: |Park& Choill2009t iPark & H wand 120091) . The study of 
substructure in clusters and of galaxy populations in different 
substructures helps us to understand the role of environmental 
effects on galaxy populations in clusters. 

To se arch for substructure in clusters we use the Mclust 
package (Fral ev & Raftervl 120061) from R, an open-source 
free statistical envir o nment developed under the GNU GPL 
dlhaka & Gentleman! 1 19961 http://www.r-project.org). 
One of our aims in this paper is to explore the possibilities of 
Mclust for the study of substructure in galaxy clusters. 

Another indicator of substructure in clusters is the deviation 
of the distribution of the peculiar velocities of galaxies in clus- 
ters (the velocities of galaxies with respect to the cluster centre) 
from Gaussian. We use several tests to analyse this distribution. 

In virialized clusters the galaxies follow the cluster poten- 
tial well. If so, we would expect that the first ranked galax- 
ies in clusters lie at the centres of groups (group haloes) and 
have l ow peculiar veloc ities (Ostriker & Tremaine 1975; Merritt 
Il984t lMalumuthlll992h . Therefore the peculiar velocity of the 
first ranked galaxies in clust ers is also an ind ication of the dy- 
namical state of the cluster dCoziol et alj|2009l) . 

Thus, in the present paper we search for subclusters in rich 
clusters of the SGW, analyse the peculiar velocities of galaxies 
in clusters, and the peculiar velocities of the first ranked galaxies. 
We study the galaxy content of substructures, using information 
about colours, luminosities, and morphological types of galax- 
ies. 

We compare our results with an N-body simulation and show 
halo merger trees to demonstrate the formation of present-day 
haloes via multiple mergers during their evolution. 

2. Data 

We use the dat a from the seventh data release of t he Sloan Digital 
Sky Survey (lAdelman-McCarthv et al.1 120081: lAbazaiian et alJ 
2009). We choose a subsample of these data from the SGW re- 
gion: 150° <R.A.< 220°, -4° < 6 < 8°, with the distance limits 
150 < D c < 300/T 1 Mpc. This region fully covers the SGW 



and also includes several poor superclusters in the foreground 
and background of the SGW, filaments which connect the SGW 
with these superclusters, and underdense regions between super- 
clusters and filaments. In this region there are 271 13 galaxies in 
the MAIN SDSS galaxy sample. 

Our next step is to determine groups of galaxies. We next de- 
scribe the sample selection and the compilation of the group cat- 
alogue; for details we refer to Tago et al. (2010). After correction 
of the galaxy magnitudes for galactic extinction we use the data 
about galaxies with the apparent r magnitudes 12.5 < r < \1 .11 . 
In addition, we have found duplicates due to repeated spec- 
troscopy for a number of galaxies in the DAS Main galaxy sam- 
ple, which had to be excluded in order to avoid false group mem- 
bers. We count as duplicate entries the galaxies which have a 
projected separation of less than 5 kpc. We excluded from our 
sample those duplicate entries with spectra of lower accuracy. 
If both duplicates had a large value of the redshift confidence 
parameter (zconf > 0.99), we excluded the fainter duplicate. 

We corrected the redshifts of galaxies for the motion rel- 
ative to the CMB and computed the co-moving distances D c 
(Martinez & Saar 2002) of galaxies using the standard cosmo- 
logical parameters: the matter density Q m = 0.27, and the dark 
energy density Q.^ - 0.73. 

We determined groups of galaxies using the friends-of- 
friends cluster analysis method introduced by Turner & Gotj 
d!976h : IZeldovich et"ai1 d!982l) : iHuchra & Gellerl dl982|) ,~aiid 
modified by us. A galaxy belongs to a group of galaxies if 
this galaxy has at least one group member galaxy closer than 
a linking length. In a flux-limited sample the density of galax- 
ies slowly decreases with distance. To take this selection effect 
properly into account when constructing a group catalogue from 
a flux-limited sample, we rescaled the linking length with dis- 
tance. As a result, the maximum sizes in the sky projection and 
velocity dispersions of our groups are similar at all distances. 
This shows that the distance-dependent selection effects have 
been properly taken into account, and that groups in our cata- 
logue form a homogeneous sample suitable for statistical stud- 
ies. However, there are exceptions: data about extremely poor 
groups with less than four member galaxies are not reliable. In 
addition, there are four exceptionally rich systems in our cat- 
alogue with more than 300 member galaxies, which all cor- 
respond to well-known nearby Abell clusters (A1656 (Coma), 
A2151/2152 (Hercules), A2197/2199 and A1367. As we see be- 
low, the clusters analysed in the present study do not belong to 
these exceptions. 

In the group catalogue the first ranked galaxy of a group is 
defined as the most luminous galaxy in the r-band. We use this 
definition also in the present paper. 

The T10 group catalogue is available in 
electronic form at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/. 

To select the galaxies and galaxy systems belonging to the 
SGW the next step was to calculate the luminosity density field 
of galaxies. 

To calculate a density field, at first we convert the spatial po- 
sitions of galaxies into spatial densities. The standard approach 
for that is to use kernel densities; we use the Bj, box spline 

B 3 (x) = — (|jc - 2| 3 - 4\x - 1| 3 + 6|jc| 3 + 4\x + 1 1 3 + \x + 2| 3 ) . 
We define the Bt, box spline kernel of a width h as 
K$\x;h)s = B 3 (x/h)(S/h), 
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Fig. 1. Distribution of groups with at least four member galaxies in the Sloan Great Wall region. Black circles correspond to the 
groups in the SGW, grey circles to the groups in the field. Circle sizes are proportional to a group's size in the sky. Crosses show the 
location of rich clusters, and numbers are their ID numbers from TableQ] 



where c5 is the grid step. This kernel differs from zero only in 
the interval x e [-2h;2h], The three-dimensional kernel func- 
K3). 

kernels: 



tion K„ ' is given by the direct product of three one-dimensional 



4 3) (r; h) s = Kf\x, h) s K^\y- h) s K^\z; h) s , 

where r = {x,y, z}. Although this is a direct product, it is 
isotropic to a good de gree. For a detailed description we refer to 
Eina sto et al.l d2007bl) . E10 and Liivamagi et al. 2010 (in prepa- 
ration). 

The ^-correction for the SDSS galaxies was calcu- 
lated with the KC ORRECT algorithm ilBlanton et all l2003at 
Blan ton & Row eis 2007). We also accepted M Q = 4.53 (in the r 
photometri c system). The evoluti on correction e was found ac- 
cording to iBlanton et al.l d2003bl) . With the velocity and radial 
dispersions from the group catalogue by T10 we supressed the 
cluster- finger redshift distortions. 

We also have to consider the luminosities of the galaxies that 
lie outside the observational window of the survey. Assuming 
that every galaxy is a visible member of a density enhancement 
(a group or cluster), we estimate the amount of unobserved lu- 
minosity and weigh each galaxy accordingly: 



L w = W L (d) L obs . 

Here WiXd) is a distance-dependent weight: 
£° L F(L)dL 



W L {d) 



rL 2 (d) 



a***™** 



(i) 



(2) 



where F(L) is the luminosity function and L\{d) and L2(d) are 
the luminosity window limits at a distance d. 

In our final flux-limited group catalogue the richness of 
groups rapidly decreases at distances D > 300/t~' Mpc due to 
selection effects. At small distances, D < 100/z -1 Mpc, the lu- 
minosity weights are large owing to the absence of very bright 
galaxies. At the distance limits used to define our sample the 
selection effects are small. 

The densities were calculated on a cartesian grid based on 
the SDSS rj, A coordinate system because it allowed the most ef- 
ficient placing of the galaxy sample cone into a box. We choose 
the kernel width h = 8/T 1 Mpc. This kernel differs from zero 
within the distance 16/z Mpc, but significantly so only in- 
side the 8/i _1 Mpc radius. We calculated the x, y, and z coor- 
dinates using the A and 77 coordinates as follows: x = -RsinA, 



y = R cos A cos 77, and z — R cos A sin 77, where R is the distance of 
the galaxy. Before extracting superclu sters we applied th e DR7 
mask assembled by Arnalte-Mur ((Martine z et al.ll2009l) to the 
density field and converted the densities into units of mean den- 
sity. The mean density is the average over all pixel values inside 
the mask. The mask is designed to follow the edges of the survey, 
and the galaxy distribution inside the mask is considered homo- 
geneous. The details of this procedure will be given in Liivamagi 
et al., 2010 (in preparation). 

Next we created a set of density contours by choosing a den- 
sity threshold and defined the connected volumes above a cer- 
tain density threshold as superclusters. Different threshold den- 
sities correspond to different supercluster catalogues. In order to 
choose proper density levels to determine the SGW and the in- 
dividual superclusters which belong to the SGW, we analysed 
the density field superclusters at a series of density levels. As a 
result we used the density level D = 4.9 to determine the indi- 
vidual superclusters in the SGW. 

The richest superclusters in the SGW are the superclusters 
SCI 126 and the supercluster SCI 111 (we use the ID num- 
bers of the superclu sters from the catalogue of superclusters by 
lEinasto et al.l (1200 ll) ). The superclust er SCI 126 a lso in cludes 
the supercluster SCI 136 from the lEinasto et al.l (12001 1) cata- 
logue. The core of the supercluster SCI 126 contains several 
rich X-ray clusters of galaxies (Belsole et al. 2004; Einasto et al. 
l2007bh . The supercluster SCI 126 is the richest in the SGW with 
a very high density core, its morphology resem bles a very rich 
and high-density multibranching filament (Eina stoet al.ll2007bl 
2008, E10). The supercluster SCI 111 is the second in richness 
in the SGW and consists of three concentrations of rich clusters 
connect ed by filaments of g alaxies (a "multispider" morphol- 
ogy, see lEinasto et al.ll2007bl) . Poor superclusters from the low- 
densi ty extension o f the S GW belong to the supercluster SCI 91 
in the lEinasto et al.! (2001 ) catalogue. These superclusters repre- 
sent different global environments in the SGW. For details about 
individual superclusters in the SGW we refer to E10. 

Then we selected rich clusters from the SGW for our study 
as follows. First, we selected all clusters of galaxies in the SGW 
with at least 75 member galaxies, altogether 7 clusters. We added 
to this sample a cluster whi ch corresponds to the Abell clus- 
ter A1773, an X-ray cluster (Bohring eTet al .1120041) . This clus- 
ter has 74 member galaxies in our catalogue. We also added 
to our analysis two clusters from the core region of the super- 
cluster SCI 126, A1650 with 48 member galaxies and A1658 
with 29 member galaxies. Both of them are X-ray clusters 
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Fig. 2. Images of some galaxies in clusters. From left to right: the first ranked galaxy in the cluster A1750 with the colour index 
g-r = 0.72, (J131236.98-011151,BRG), a red spiral galaxy in the cluster A1066 (J10351. 33+051 130.8, g - r = 0.75, BRG), two 
galaxies in the cluster A1658 (J1301 16.68-032823.6 with g - r = 0.67, BRG, and J1301 16.09-032829.1 with g - r = 0.54), and a 
red spiral galaxy in the cluster A1663 (J 130 152.6-0240 12.3 with g - r = 0.75, BRG). 

Table 1. Rich clusters in the SGW. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


SclID 


No 


ID 

1 


Dist. 
h~ l Mpc] 


R.A. 
[deg] 


Dec. 

[deg] 


L tot 
10 10 /T 2 L o 


0"V 

kms~ l 


N gal 


N BRG 


N BRS 


SCI 91 


1 


1066 


206.6 


159.8 


5.2 


0.854E+02 


758 


78 


15 


8 




2 


1205 


228.4 


168.6 


2.3 


0.108E+03 


624 


105 


11 


1 


SCI 111 


3 


1424 


226.0 


179.3 


5.1 


0.968E+02 


656 


78 


12 


4 




4 


1516 


230.3 


184.4 


3.7 


0.119E+03 


834 


104 


15 


3 


SCI 126 


5 


1650x 


249.8 


194.7 


-1.7 


0.706E+02 


466 


48 


8 


3 




6 


1658x 


253.4 


195.3 


-3.5 


0.445E+02 


441 


29 


4 


1 




7 


1663x 


246.7 


195.6 


-2.6 


0.104E+03 


530 


78 


12 


2 




8 


1750x 


256.4 


202.7 


-1.9 


0.190E+03 


660 


117 


17 


2 




9 


1773x 


230.6 


205.6 


2.3 


0.918E+02 


646 


74 


11 


4 




10 


1809 


236.5 


208.3 


5.2 


0.118E+03 


649 


90 


12 


3 



Columns in the Table are as follows: Cols. 1-3: Supercluster ID, cluster No and Abell ID (x denotes X-ray clusters), Cols. 4- 



R.A. and Dec. of the cluster, Col. 7: total luminosity of the supercluster L[ (, Col. 8: cry - rms velocity, Cols. 9-11: 
^BRG ~~ numrjer of bright red galaxies, Ngj^g - number of spirals among the bright red galaxies. 



N 



gal 



6: comoving distance, 
- number of galaxies, 



(iBohringer et al.ll2004l) . The reason is that we wanted to study in 
detail all rich clusters in the core region of this supercluster. Our 
earlier analysis (E10) has shown that there are several differences 
in the properties of groups and of the galaxy content between the 
core of the supercluster SCI 126 and other superclusters in the 
SGW: groups in the core of the supercluster SCI 126 are richer 
than in other superclusters of the SGW, and the fraction of red 
galaxies there is larger than in the outskirts of this superclus- 
ter or in other superclusters in the SGW. Now we compared the 
properties of rich clusters in the core of the supercluster SCI 126 
with the properties of rich clusters from other superclusters in 
the SGW to see whether the dynamical state of clusters differs. 

Thus, our sample of rich clusters consists of ten clusters and 
includes all very rich clusters in the SGW, all X-ray clusters and 
all rich clusters from the core of the supercluster SCI 126. All 
these clusters correspond to Abell clusters; their Abell ID and 
other data are given in Table [TJ We show the sky distribution of 
groups in the SGW in Fig. [TJ In this figure clusters are marked 
with numbers from Table QJ The clusters 1 and 2 are located 
in the supercluster SCI 91, the clusters 3 and 4 in the super- 
cluster SCI 111, the clusters 5-8 in the core of the supercluster 
SCI 126 and the clusters 9 and 10 in the outskirts of the super- 
cluster SCI 126. 

Our clusters are chosen from a narrow distance interval 
(Table [TJ, here our sample is almost volume-limited. Absolute 
magnitude limits of galaxies in groups in the superclusters 
SCI 91 and SCI 111 are approximately -19.20, and in the su- 
percluster SCI 126 -19.40. 



3. Methods 

3. 1 . Substructure in rich clusters 

We studi ed the substructure o f clusters with the R package 
Mc/«sf( Fr alev & Raftervll2006l) . This package searches for an 
optimal model for the clustering of the data among the models 
with varying shape, orientation, and volume, and for the optimal 
number of components or substructures, using multidimensional 
normal mixture modeling. Below we will use Mclust to calcu- 
late the number of components in clusters using the data about 
the sky position and peculiar velocities of the cluster member 
galaxies. To see how these components are populated with galax- 
ies of different colour and luminosity we applied Mclust to an 
extended dataset including the colours and luminosities of the 
cluster members. Mclust calculates for every galaxy the proba- 
bilities to belong to any of the components. The uncertainty of 
classification is defined as 1. minus the highest probability of a 
galaxy to belong to a component. The mean uncertainty for the 
sample is used as a statistical estimate of the reliability of the 
results. 

When we searched for substructure in 3D, we assumed that 
the line-of-sight positions of galaxies are given by their pecu- 
liar velocities, because we had no other hypothesis that could be 
used. However, while peculiar velocities might carry some dis- 
tance information, the only fact we can safely assume is that the 
cluster galaxies are located inside the cluster volume. This is the 
reason why we studied the velocity distribution separately. We 
also tested how the possible errors in the line-of-sight positions 
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Table 2. Dynamical properties of rich clusters. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


SclID 


No 


ID 


Ncomp 


v pec 


Vrpec 


PSW 


k 


Pk 


sk 


Psk 


SCI 91 


1 


1066 


1 


-312 


0.41 


0.074 


2.34 


0.139 


-0.243 


0.528 




2 


1205 


5 


49 


0.08 


0.247 


2.37 


0.099 


-0.069 


0.841 


SCI 111 


3 


1424 


2 


-825 


1.26 


0.173 


2.51 


0.387 


0.383 


0.339 




4 


1516 


1 


-58 


0.07 


0.535 


2.84 


0.948 


-0.292 


0.403 


SCI 126 


5 


1650x 


1 


-189 


0.41 


0.110 


2.80 


0.946 


0.563 


0.263 




6 


1658x 


1 


-463 


1.05 


0.937 


2.37 


0.583 


0.048 


0.936 




7 


1663x 


2 


-841 


1.59 


0.372 


2.57 


0.493 


-0.228 


0.562 




8 


1750x 


5 


-724 


1.10 


0.115 


2.51 


0.243 


-0.265 


0.422 




9 


1773x 


1 


-76 


0.12 


0.572 


3.24 


0.436 


-0.083 


0.835 




10 


1809 


2 


-115 


0.18 


0.802 


2.73 


0.780 


0.053 


0.887 



Columns in the Table are as follows: 1-3: Supercluster ID, cluster No and Abell ID, Col. 4: number of components, Col. 5: peculiar velocity of 
the first ranked galaxy (bis -1 ), Col. 6: normalised peculiar velocity of the first ranked galaxy, Vpec/cv', Col. 7: p-value of the Shapiro-Wilk test, 
Cols. 8 and 9: kurtosis k and the p-value pj, for kurtosis for the peculiar velocities, Cols. 10 and 11: skewness sk and the p- value p s j, for skewness 
for the peculiar velocities. 

Table 3. Galaxy populations in rich clusters. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


SclID 


No 


ID 


(g ~ 'Omean 


F red 


F ell 


Frs 


slope re( j 


SD red 


SD Ered 


SD Sred 


SCI 91 


1 


1066 


0.73 


0.73 


0.36 


0.58 


-0.024 


0.043 


0.024 


0.053 




2 


1205 


0.68 


0.62 


0.38 


0.45 


-0.012 


0.030 


0.027 


0.033 


SCI 111 


3 


1424 


0.73 


0.71 


0.44 


0.49 


-0.007 


0.046 


0.033 


0.057 




4 


1516 


0.76 


0.85 


0.56 


0.34 


-0.018 


0.032 


0.032 


0.033 


SCI 126 


5 


1650x 


0.78 


0.94 


0.56 


0.40 


-0.010 


0.039 


0.036 


0.042 




6 


1658x 


0.74 


0.83 


0.58 


0.29 


-0.026 


0.029 


0.031 


0.024 




7 


1663x 


0.75 


0.83 


0.49 


0.43 


-0.007 


0.037 


0.032 


0.042 




8 


1750x 


0.78 


0.91 


0.62 


0.33 


-0.013 


0.042 


0.041 


0.046 




9 


1773x 


0.73 


0.77 


0.53 


0.41 


-0.011 


0.040 


0.032 


0.048 




10 


1809 


0.75 


0.84 


0.51 


0.43 


-0.015 


0.029 


0.027 


0.032 



Columns in the Table are as follows: Cols. 1-3: Supercluster ID, cluster No and Abell ID, Col. 4: (g - r) m ean - the mean g — r colour of galaxies, 
Cols. 5-11: fraction of red galaxies, fraction of elliptical galaxies, fraction of spirals among red galaxies, and the slope and the rms scatter of the 
red sequence in the colour - magnitude diagram (the scatter of all red galaxies, and that of red ellipticals and red spirals). 



of galaxies affect the results of Mclust randomly shuffling the pe- 
culiar velocities of galaxies 1000 times and searching each time 
for the substructure with Mclust. The number of the components 
found by Mclust remained unchanged, which demonstrated that 
the results of Mclust are robust. 

To analyse the distribution of the peculiar velocities of galax- 
ies in clusters we used several ID tests. We tested the hy- 
pothesis about the Gaussian distribution of the peculiar veloc- 
ities of galaxie s in clusters with the Shapiro-Wilk normality test 
(Shapiro 1965). We used this test because it is considered the 
best for small samples. We also calculated the kurtosis and the 
skewness of the peculiar velocity distributions, and used these 
to test for the normality of the distributions with the Anscombe- 
Glynn test for the kurtosis ( Anscombe & Glvnn 1983) and with 
the D'Agostino test for the skewness dD'Agost ino 1970), us- 
ing the R package moments by L. Komsta and F. Novomestky 
(http : //www . r-pro j ect . org). 

The distribution of the peculiar velocities of galaxies in clus- 
ters can be quantified wi t h the usual normal mixture models 
(e.g. | Ashman et alJ (11994 : iBoschin et al.1 (120061) : iBarrena et alJ 
(2007). In the R environment there is a s uitable packag e to study 
these mixture models, called flexmix dLeischl 120041) . Flexmix 
fits a user-specified number of Gaussians to the velocity data. 



However, subclusters in clusters sometimes do not differ in pe- 
culiar velocities, and in some clusters the tails of the velocity 
distribution are not well determined by flexmix. Thus we fitted 
the distributions of the peculiar velocities with normal mixture 
models ourselves, taking into account the substructure informa- 
tion provided by Mclust and by the 3D distribution of galaxies 
that we will show below. 

The data characterizing the clusters' dynamical state are 
given in Table|2] where we present the number of components in 
clusters determined by Mclust, the peculiar velocities of the first 
ranked galaxies, the results of the Shapiro-Wilk test, the kurtosis 
and skewness, and their normality test /^-values for the peculiar 
velocity distributions. 

3.2. Galaxy populations in rich clusters 

To analyse the galaxy content of clusters, we used the g—r colour 
of galaxies and their absolute magnitude in the r-band M r . We 
divided the galaxies by colour into the red and blue populations 
using the colour limit g — r — 0.7 (red galaxies have g - r > 0.7) 
and excluded from the analysis the galaxies with g - r > 0.95 
(such a high value indicates possible errors in photometry, often 
due to overlapping images of galaxies). This limit depends on the 
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luminosity of galaxies; because there are no very faint galaxies 
in our sample, we will use this simple approach. We calculated 
the fraction of red galaxies in clusters, the mean colour of the 
cluster, and the rms scatter and slope of the red sequence in the 
cluster colour-magnitude diagram (for all galaxies, for red ellip- 
ticals, and for red spirals). 

We also studied the population of bright red galaxies (BRGs, 
the galaxies that have the GALAXY_RED flag in the SDSS 
database) in the clusters. The BRGs are nearby (z < 0-15 , cut I 
LRGs (Eisens tein et al.l l2001). Because [Eisenstein et al.l (12001 
warn that the sample of nearby LRGs may be contaminated by 
galaxies of lower luminosity, we choose to call them BRGs. The 
BRGs are similar to the LRGs at higher redshifts. Nearby bright 
red gala xies do not form an ap proximately volume-limited pop- 
ulation (Eisenstein et al. 2001) but they are yet the most bright 
and the most red galaxies in the SGW region (see also E10). 

We determined the morphological type of galaxies in clusters 
with the SDSS v isual tools. Follow ing the classification used in 
the ZOO project dLintott et al.ll2.Q08h . we denote elliptical galax- 
ies as of type 1, spiral galaxies of type 4, those galaxies which 
show a signs of merging or other disturbancies as of type 6. 
Types 2 and 3 mark clockwise and anti-clockwise spiral galaxies 
with visible spiral arms. Fainter galaxies, for which the morpho- 
logical type is difficult to determine, are of type 5. We calculate 
the fraction of elliptical galaxies and of red spirals in clusters. A 
word of caution is needed: faint galaxies are sometimes difficult 
to classify reliably. Some red galaxies classified as spirals may 
actually be SO galaxies. We present images of some galaxies in 
clusters in Fig. [2] All these galaxies are classified as BRGs (an 
exception is a blue galaxy in the third image from the left, which 
may be interacting with a companion galaxy). 

The data for cluster galaxy populations are given in Table [3] 
Images of all galaxies in our clusters are shown on the web page 
http : //www . aai . ee/~maret/SGWcl . html, 

In the following analysis of individual rich clusters we 
present for each cluster the sky distribution of galaxies along 
with the density contours, the 3D vizualisation plot of galaxies 
(to generate these figures, we used the R package scatterplot3D 
by Ligges & Machler (2003)), the histograms of the peculiar ve- 
locities for galaxies of a different morphological type and for the 
BRGs, and the colour-magnitude diagrams. We have added the 
Gaussians for the components to the histograms of the peculiar 
velocities of all galaxies. For some clusters we include additional 
figures to show their substructure. 

In the figures of the sky distribution of galaxies in clusters the 
2D density contours are calculated with Mclust for clusters with 
more than one component, and with the ^package Kernsmooth 
for clusters with one component only. 



4. Results 

4.1. Clusters in the superclusters SC1 111 and SCI 91 

4.1.1. The cluster A1 066 

The first rich cluster we study is the cluster A 1066, which is 
located in the supercluster SCI 91. Figure [3] (left panel) shows 
that the 2D density contours of this cluster are quite smooth el- 
lipsoids. According to Mclust, this cluster consists of one com- 
ponent, the mean uncertainty of the classification is 1.89 • 10~ 2 . 
However, the 3D distribution of galaxies in this cluster Fig. [3] 
(right panel) is complex - the shape resembles an hourglass, 
where some galaxies with negative peculiar velocities form a 
separate component. 
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Fig. 3. 2D and 3D view of the cluster A 1066 (SCI 91). Open cir- 
cles correspond to elliptical galaxies, crosses - to spiral galaxies, 
filled circles to BRGs, and a star to the first ranked galaxy in the 
cluster. 
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Fig. 4. Left panel shows R.A. vs. Dec, the right panel RA. vs. 
peculiar velocities of galaxies (in kms~ l ) in the cluster A1066. 
Different symbols correspond to three different velocity com- 
ponents (see text), ellipses show the covariances of the compo- 
nents. 
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Fig. 5. Distribution of the peculiar velocities of galaxies in the 
cluster A1066. Top down: all galaxies, elliptical galaxies, spiral 
galaxies (the grey histogram shows the red spirals), and BRGs. P 
denotes the probability histogram, F the galaxy number counts. 
The small dash in the upper panel indicates the peculiar velocity 
of the first ranked galaxy. 



Figure |U a so-called classification plot produced by Mclust, 
shows three components in the distribution of the peculiar ve- 
locities of galaxies. The peculiar velocities with negative val- 
ues correspond to the lower part of the "hourglass". The left 
panel of this figure shows that the different components in ve- 
locity are not separated in the sky distribution, which is why 
Mclust determines one component only for this cluster as given 
in Table [2] We see three components in the distribution of the 
peculiar velocities of galaxies in this cluster, and the p-value of 
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Fig. 6. Colour-magnitude diagram of galaxies in the cluster 
A1066. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 



the Shapiro- Wilk test for this cluster is 0.074, the lowest among 
our sample (the lower the /rvalue, the higher is the difference to 
the Gaussian null hypothesis). 

The histograms of the peculiar velocities for all galaxies and 
separately for the elliptical and spiral (and red spiral) galax- 
ies and for the BRGs in A1066 are shown in Fig. [5] At the 
peculiar velocities v pec > lOOOkm/s we see a small separate 
component of galaxies consisting mostly of spirals (including 
red spirals), although there are also some elliptical galaxies and 
one BRG among them; this system is also seen in Fig. [4] In 
the main component of this cluster the galaxies have peculiar 
velocities 1000 > v pec > -500 km/s, the first ranked galaxy 
of the cluster also belong here. At negative peculiar velocities 
(Vp e c < -500km/s) there is another extension (a lower part of the 
"hourglass" in Fig. [3] right panel) of galaxies, consisting again 
mostly of spirals, including red spirals and five BRGs. Eight 
of sixteen BRGs in this cluster have been classified as spirals 
(Table [3]), five of them are located in this extension. The distri- 
bution of the peculiar velocities of elliptical galaxies shows a 
concentration towards the cluster centre, while the distribution 
of the peculiar velocities of spiral galaxies is smoother. 

In the colour-magnitude diagram for A 1066 (Fig. [6]) we show 
with separate symbols elliptical and spiral galaxies and BRGs. 
The slope and the rms scatter of the red sequence in the colour- 
magnitude diagram are presented in Table[3] 73% of the galaxies 
in this cluster are red, more than half of the red galaxies are spi- 
rals (or SOs). The colour-magnitude diagram shows that the scat- 
ter of colours of red spirals (both the red spiral BRGs and fainter 
red spirals) is larger than the scatter of colours of red ellipti- 
cal galaxies (two faint blue galaxies in this cluster are classified 
as blue ellipticals), which suggests that red elliptical galaxies in 
this cluster form a more homogeneous population than red spiral 
galaxies. High negative and high positive peculiar velocities of 
some spiral galaxies and some BRGs in this cluster suggest that 
they dynamically form separate components of the cluster. 

Summarizing, this cluster, classified by Mclust as a one- 
component cluster with a rather smooth sky distribution of 
galaxies has indeed a substructure, which suggests that this clus- 
ter has not yet completed its formation. 



4.1.2. The cluster A1 205 

The sky distribution of galaxies in the cluster A 1205 in the 
supercluster SCI 91 (Fig. |7J left panel) shows three different 
components. According to Mclust, this cluster even consists of 
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Fig.7.2D and3D views of the cluster A 1205 (SCI 91). Symbols 
are the same as in Fig. [3] 
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Fig. 8. Distribution of the peculiar velocities of galaxies in the 
cluster A 1205. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. The small dash in the upper panel indicates 
the peculiar velocity of the first ranked galaxy. 
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Fig. 9. Left panel shows R.A. vs. Dec, the right panel the R.A. 
vs. the peculiar velocities of galaxies (in kms~ l ) in the cluster 
A 1205. Different symbols correspond to different velocity com- 
ponents. 



five components, the mean uncertainty of the classification is 
1.62 • 10~ 3 . The reason for that is the distribution of the pecu- 
liar velocities of galaxies in each component that is shown in 
Fig. [9] This figure shows that the galaxies in a left component 
are divided into three according to their peculiar velocities; part 
of them have high positive peculiar velocities and others mostly 
negative peculiar velocities. There is also a small subcluster of 
galaxies with very low peculiar velocities. 

In another big component the galaxies have mostly positive 
peculiar velocities, in the third one - negative peculiar velocities 
(Fig. [8]). The shape of the 3D distribution of galaxies Fig.|7](right 
panel) resembles an hourglass, this time seen at an angle. 

Even more interesting is the structure of this cluster if we 
plot the sky distribution with the density contours for red galax- 
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Fig. 10. Distribution of red galaxies (filled circles) in the cluster 
A 1205. The star denotes the first ranked galaxy. 
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Fig. 11. Colour-magnitude diagram of galaxies in the cluster 
A 1205. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 



ies only (Fig.fTOt. This figure shows that one component in this 
cluster mostly consists of red galaxies (there are both elliptical 
and spiral galaxies among them, compare Fig. UJ left panel and 
Fig.ES. 

The first ranked galaxy in this cluster is located almost at the 
centre of the sky distribution and has a very low peculiar velocity 
(the lowest among the clusters in our sample, see Table|2|. 

The histograms of the peculiar velocities of galaxies in 
A1205 in Fig. [8] show a large fraction of red spirals with neg- 
ative peculiar velocites, which belong to the component with red 
galaxies. Some BRGs are also located there; in addition, we see 
several BRGs in other components, even at the edges of the clus- 
ter (Fig. |7]). The peculiar velocities of galaxies from the three 
different components seen in the sky distribution partly overlap, 
so these components are not very clearly seen in the histograms. 

In the colour-magnitude diagram for the A1205 (Fig.fTTI) we 
see a bimodal distribution of galaxies, where the red and blue 
sequences are both clearly seen. This cluster has the largest frac- 
tion of blue galaxies in our sample (32%). The colour-magnitude 
diagram shows that the scatter of colours of red spirals is larger 
than the scatter of colours of red elliptical galaxies (the rms scat- 
ters are 0.033 and 0.027, correspondingly, for red spirals and 
for red ellipticals), although the difference is not as pronounced 
as in the cluster A 1066. In this cluster there are four faint blue 
galaxies classified as blue ellipticals. 

The presence of substructures in this cluster, which are popu- 
lated by galaxies of different properties, suggests that the cluster 
A 1205 consists of at least three merging components. 
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Fig. 12. 2D and 3D views of the cluster A1424 (SCI 111). 
Symbols are the same as in Fig. [3] 
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Fig. 13. Sky distribution of elliptical (left) and spiral (right) 
galaxies in the cluster A1424. The 2D density contours are cal- 
culated using the data about elliptical or spiral galaxies, respec- 
tively. Filled circles in the left panel denote elliptical galaxies, in 
the right panel spiral galaxies. The star denotes the first ranked 
galaxy. 



4.1.3. The cluster A1 424 

The next cluster under study, the cluster A1424, is located in the 
second rich supercluster in the SGW, SCI 111. The sky distri- 
bution of galaxies in this cluster (Fig. [12] left panel) shows two 
components also confirmed by Mclust. The mean uncertainty of 
the classification of galaxies in the cluster A1424 is 2.92 • 10~ 4 . 

The 3D distribution of galaxies (Fig. [12] right panel) shows 
that in one of these components galaxies have mostly positive 
peculiar velosities. Figure [T3] shows that according to the sky 
distribution of galaxies in the cluster A1424, the components 
as delineated by elliptical and spiral galaxies are different - the 
two components are more clearly separated in the distribution of 
spiral galaxies. 

The distribution of the peculiar velocities of galaxies in the 
cluster A1424 (Fig. RBI has a minimum near the clusters cen- 
tre, which is seen in the distribution of the peculiar velocities 
of both elliptical and spiral galaxies. The elliptical galaxies are 
more concentrated towards the centre of the cluster than the spi- 
ral galaxies, which is why the distribution of the peculiar veloci- 
ties is smoother. In the same time, the distribution of the peculiar 
velocities of red spirals also shows two separate components. 
The BRGs are mostly located in the component with smaller 
right ascensions, several of which are spirals. The first ranked 
galaxy of this cluster is located almost in the centre of this com- 
ponent. 

In this cluster 73% of all galaxies are red. In the colour- 
magnitude diagram (Fig.fTBTl the scatter of red spiral galaxies is 
about two times larger than the scatter of red elliptical galaxies 
(the rms scatters are 0.033 and 0.057, correspondingly, for red 
ellipticals and for red spirals). This is caused by galaxies in the 
outer parts of the cluster; it is possible that these galaxies have 
only recently joined the cluster. 
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Fig. 14. Distribution of the peculiar velocities of galaxies in the 
cluster A1424. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. The small dash in the upper panel indicates 
the peculiar velocity of the first ranked galaxy. 
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Fig. 15. Colour-magnitude diagram for galaxies in the cluster 
A1424. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 
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Fig. 17. Left panel shows R. A. vs. Dec, the right panel the R A. 
vs. the peculiar velocities of galaxies (in kms~ l ) in the cluster 
A1516. Different symbols correspond to different velocity com- 
ponents. 
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Fig. 18. Distribution of the peculiar velocities of galaxies in the 
cluster A1516. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. The small dash in the upper panel indicates 
the peculiar velocity of the first ranked galaxy. 




Fig. 16. 2D and 3D views of the cluster A1516 (SCI 111). 
Symbols are the same as in Fig. [3] 



4.1.4. The cluster A1 51 6 

The cluster A1516 is located in the high-density core of the su- 
percluster SCI 111. The sky distribution of galaxies in this clus- 
ter shows very regular contours (Fig. [16J left panel). According 
to Mclust, this cluster has one component, the mean uncertainty 
of the classification of galaxies is 7.71 ■ 10~ 4 . 

The 3D distribution of galaxies shows that in the sky galaxies 
with negative peculiar velocities are mostly located in the central 
part of the cluster (Fig. [T6] right panel, and Fig. [T7b . where they 



form an extension with the outer parts delineated mostly by spi- 
ral galaxies. This component is seen also in Fig. [18] According 
to the peculiar velocities, the elliptical galaxies in this compo- 
nent are concentrated towards the centre of the cluster. There 
is another peak in the distribution of the peculiar velocities of 
elliptical galaxies closer to the cluster centre; a signature of a 
past merger? The distribution of the peculiar velocities of spiral 
galaxies is less peaky. The distribution of the peculiar velocities 
of red spiral galaxies shows a concentration at the centre, at low 
peculiar velocites. There is also a small subcluster of galaxies 
with high positive peculiar velocities. The gradient of the pecu- 
liar velocities of galaxies in A1516 hints that this cluster may be 
rotating. With one exception, BRGs populate central parts of the 
cluster, and the first ranked galaxy of the cluster has a very low 
peculiar velocity, v pec - 58 km/s. The /?-value of the Shapiro- 
Wilk test is 0.535, which confirms the Gaussian distribution of 
the peculiar velocities of galaxies in this cluster. 

The colour-magnitude diagram in Fig. [19] shows that the 
scatter of colours of red ellipticals and red spirals is almost equal. 
The fraction of red galaxies in this cluster is very high (85%, 
Table [3J, only 3 of 15 BRGs in this cluster are classified as spi- 
rals. 
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Fig. 19. Colour-magnitude diagram of galaxies in the cluster 
A1516. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 
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Fig. 20. 2D and 3D views of the cluster A1650. Symbols are the 
same as in Fig. [3] 



4.2. Clusters in the core of the supercluster SCI 126 

4.2.1. The cluster A1 650 

The first cluster we study in the core of t he supercluster SCI 126 
is the cluster A165 0, an X-ray cluster ([Bohringer et al. 2004; 
IUdomprasertl2 004). This cluster is the second-poorest cluster in 
our sample with only 48 member galaxies in our catalogue. The 
sky distribution of galaxies in this cluster (Fig. [20] left panel) 
shows two concentrations. The 3D view and the velocity his- 
tograms of galaxies of different type (Fig. [20] right panel, and 
Fig. l2"TT i show an almost separate component of galaxies with 
positive peculiar velocities, all but one galaxies in this com- 
ponent are elliptical. In the sky distribution this component is 
located almost in the centre of the cluster, thus Mclust, using 
the data about all the galaxies, finds that this cluster consists of 
one component, with the mean uncertainty of the classification 
0.008. Red spiral galaxies and BRGs ( with one exception) popu- 
late the main component of the cluster. lUdompra sert ( 2004J) de- 
scribes this cluster as a compact X-ray source, possibly located 
at a cold spot in the CMB. 

However, the sky distribution of only elliptical or only spiral 
galaxies (Fig.l22l clearly shows two components in the distribu- 
tion of elliptical and spiral galaxies, in both cases in a different 
manner. Thus we approximated the distribution of the peculiar 
velocities of galaxies with three Gaussians. These are better seen 
in the distribution of the peculiar velocities of elliptical galax- 
ies. The distribution of the peculiar velocities of spiral galax- 
ies is different. There is one spiral galaxy only in one compo- 
nent, which has the highest peculiar velocity in this component. 
Mclust confirms the components in the distribution of elliptical 
and spiral galaxies. The first ranked galaxy of this cluster is lo- 
cated at the edge of the elliptical galaxy component, but in the 
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Fig. 21. Distribution of the peculiar velocities of galaxies in the 
cluster A 1650. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. The small dash in the upper panel indicates 
the peculiar velocity of the first ranked galaxy. 
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Fig. 22. Distribution of elliptical (left panel) and spiral (right 
panel) galaxies in the cluster A 1650. Filled circles in left panel 
denote elliptical galaxies, in right panel they show spiral galax- 
ies. The star is the first ranked galaxy. 



centre of one spiral galaxy component; the peculiar velocity of 
the first ranked galaxy v pec — -189 km/s. 

The fraction of red galaxies in this cluster is very high, some 
spiral galaxies are also red. Owing to the small number of galax- 
ies in this cluster the high fraction of red galaxies is only ap- 
proximate. The first ranked galaxy in this cluster is also a red 
spiral galaxy, as seen in the colour-magnitude diagram (Fig.l23l. 
In this diagram the rms scatter of colours of the red sequence is 
0.036 for red ellipticals and 0.042 for red spirals. 

Our analysis hints at a possibility that this cluster consists 
of two main components, some galaxies belong to a third sub- 
group with high positive peculiar velocities. Because of the small 
number of galaxies in this cluster our results about substructures 
should be taken as a suggestion only. 

4.2.2. The cluster A1 658 

The next cluster we study in the core of the supercluster SCI 126 
is the cluster A1658, the poorest cluster in our sample with only 
29 member galaxies. The sky distribution of galaxies in this clus- 
ter (Fig. [24] left panel) shows three irregularly located concen- 
trations. According to Mclust, using the data about all the galax- 
ies, this cluster has one component, the mean uncertainty of the 
classification of galaxies is less than 10~ 5 . 

The distribution of the peculiar velocities of galaxies in 
A1658 is shown in Fig. [25] The distribution of the peculiar ve- 
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Fig. 23. Colour-magnitude diagram of galaxies in the cluster 
A 1650. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 




Fig. 24. 2D and 3D view of the cluster A1658. Symbols are the 
same as in Fig. [3] 



locities of spiral galaxies shows two separate components, while 
the distribution of the peculiar velocities of elliptical galaxies has 
a central concentration, and two elliptical galaxies have high pe- 
culiar velocities. These distributions may be a signature of merg- 
ing, which is more clearly seen in the distribution of the peculiar 
velocities of spiral galaxies. The p-value of the Shapiro-Wilk test 
is p — 0.937, the highest among our cluster sample, confirms the 
Gaussianity of the velocity distribution. The peculiar velocity of 
the first ranked galaxy in this cluster v pec - -463 km/s; this 
galaxy is located in the richest component of this cluster. There 
are BRGs in both components, and most of the red spiral galax- 
ies belong to one component. 

In the colour-magnitude diagram the scatter of colours of red 
galaxies in this cluster is very small, the fraction of red galaxies 
is high. Surprisingly, the scatter of colours of red spirals is even 
smaller than that of red ellipticals (with the rms values of 0.024 
and 0.031, respectively). If the bimodal distribution of the pecu- 
liar velocities of galaxies in this cluster is an indication about a 
recent merger of two clusters into one, then the colour evolution 
of galaxies was probably already finished in these clusters before 
merging. Again, as we also mentioned with regard to the cluster 
A 1650, the number of galaxies in this cluster is very small and 
our results concerning this cluster should be taken as suggestion. 
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Fig. 25. Distribution of the peculiar velocities of galaxies in the 
cluster A1658. Top down: all galaxies, elliptical galaxies, spi- 
ral galaxies (grey - red spirals), and BRGs. The small dash in 
the upper panel indicates the peculiar velocity of the first ranked 
galaxy. 
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Fig. 26. Colour-magnitude diagram of galaxies in the cluster 
A1658. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 
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Fig. 27. 2D and 3D views of the cluster A 1663. Symbols are the 
same as in Fig. [3] 



4.2.3. The cluster A1 663 

Figure [27] shows that the distribution of galaxies in the cluster 
A 1663 has a concentration where mostly galaxies with positive 
peculiar velocities are located. These galaxies form a separate 
component in the velocity distribution (Fig. [29l . According to 
Mclust, this cluster has two components with the mean uncer- 
tainty of the classification 9.54 • 10~ 3 . 

The distribution of the peculiar velocities of all galaxies in 
A1663 in Fig. [29] shows three components. A concentration of 



the peculiar velocities near the cluster centre is mostly due to 
elliptical galaxies. The distribution of the peculiar velocities of 
spiral galaxies, especially that of red spirals, shows three sub- 
systems, one of them corresponding to the component with high 
negative peculiar velocities. The p-value of the Shapiro- Wilk 
test is p — 0.37, so it does not confirm the non-Gaussianity of 
the velocity distribution. The BRGs in this cluster have peculiar 
velocities v pec < 500 km/s, and in the sky distribution they are 
located at higher densities (according to the density contours in 
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Fig. 29. Distribution of the peculiar velocities of galaxies in the 
cluster A1663. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. The small dash in the upper panel indicates 
the peculiar velocity of the first ranked galaxy. 
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Fig. 30. Colour-magnitude diagram of galaxies in the cluster 
A 1663. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 



Fig. [27] where the first ranked galaxy is also located. The first 
ranked galaxy in this cluster has a high negative peculiar veloc- 
ity (vpec — -841 km/s), sugg esting that this cluster is not viri- 
alized yet. iBurgett et al.l (120041) suggest from the data from the 
2dFGRS that the velocity gradient in this cluster is an indication 
that the cluster may be rotating. 

The fraction of red galaxies in the cluster A 1663 is high, 
83%, and 43% of red galaxies are spirals. In the colour- 
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Fig. 31. 2D and 3D views of the cluster A1750. Symbols are the 
same as in Fig. [3] 



magnitude diagram the red and blue galaxies in this cluster are 
clearly separated. The scatter of colours of red galaxies in the 
cluster A1663 is larger than in the cluster A1658, with the rms 
value 0.037, and the scatter of colours of red ellipticals is larger 
than that of red spirals (with the rms values of 0.032 and 0.042, 
correspondingly ) . 

4.2.4. The cluster A1 750 

The next cluster, which is the richest cluster in our sample with 
117 member galaxies in our group catalogue, is the well-known 
binary merging Abell cluster A1750. 

Mclust determined even five components in this cluster. In 
Fig.[5T|(left panel) we mark them with numbers in an increasing 
order of the right ascencion. The mean uncertainty of the classi- 
fication of galaxies in the cluster A1750 is 9.8 • 10~ 3 . 

The 3D distribution of galaxies in this cluster (Fig.[3T] right 
panel) shows two almost separate components. In one of them 
most galaxies have negative peculiar velocities with a mean 
value v pec - -600 km/s. This component, marked as compo- 
nent 4 in the left panel of Fig. [31] is well seen in the classifi- 
cation diagram generated with Mclust, Fig. [32] The first ranked 
galaxy of the cluster A 1750 is also located in this component 
and has the peculiar velocity v pec = -724 km/s. In Fig. [2] (left 
panel) we show an image of the first ranked galaxy in this clus- 
ter. This galaxy, a BRG, is embedded in a common red halo with 
its nearest companion galaxy which is the second in brightness 
in this cluster. These galaxies may form a merging system. The 
histograms of the peculiar velocites (Fig.[33l show that this com- 
ponent is better seen in the distribution of spiral galaxies than in 
the distribution of elliptical galaxies. Some BRGs are also lo- 
cated in component 4. 

In Fig. [5T|the 2D density contours show that the sky den- 
sity of galaxies in the cluster A 1750 is the highest in component 
3. The mean peculiar velocity of galaxies in this component is 
Vpec = -72 km/s. Several BRGs are located here. 

In another component, component 5, galaxies have positive 
peculiar velocities, v pec = 612 km/s. There are two BRGs in this 
component, and some elliptical and spiral galaxies. 

In two components, 1 and 2, the mean peculiar velocity of 
galaxies is almost equal, v pec ~ 220 km/s, and the 2D density of 
galaxies decreases in these components as we move farther away 
from component 3 with a high 2D galaxy density. Figure [3T| 
shows that also some BRGs are located in these components, 
in the outer parts of the cluster. 

The distribution of the peculiar velocites of red spiral galax- 
ies (Fig. [33} follows that of all spirals (of course, most spirals in 
this cluster are red). 

The fraction of red galaxies in this cluster is very high (9 1 %), 
and the rms scatter of colours of red elliptical and red spiral 
galaxies is almost equal, about 0.04. The colours of galaxies in 
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Fig. 32. Left panel shows R.A. vs. Dec, the right panel the 
R.A. vs. peculiar velocities of galaxies (in kms~ l ) in the cluster 
A 1750. Different symbols correspond to different velocity com- 
ponents. 




-2000 -1000 0.0 
, n v(pec)(km/s) 



_h: 



LB_ 



m 



-2000 

10 



-1000 0.0 

v(pec)(km/s) 



Jib 



0.0 

v(pec)(km/s) 



Fig. 33. Distribution of the peculiar velocities of galaxies in the 
cluster A1750. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. 
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Fig. 34. Colour-magnitude diagram of galaxies in the cluster 
A 1750. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 



different components are rather similar. In this respect the clus- 
ter A1750 differs from the cluster A1205, where one component 
consists mostly of red galaxies. 

The structure of the cluster A 1750 has been analysed in sev- 
eral o t her studies using op t ical and X-ray data (Donnelly et al. 
2001; iBelsole et all 12004 iBurgett et a l. 2004; Hwang &Leel 



2009). The comparison of substructures shows that our compo - 



nent 3 coincides wit h the component C bv lHwang & Led (12009). 
IBelsole et al.l (|2004) finds an X-ray peak at this location and sug- 
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Fig. 35. 2D and 3D views of the cluster A 1773. Symbols are the 
same as in Fig. [3] 



gest that this component suffered a merger or interaction in the 
past 1-2 Gyrs. Our component 4 is the same as the one deter- 
mined by iHwang & Led d2009l) as component N, who suggest 
that the components 3 and 4 have started interacting. Component 
4 coin cides with another X-ray peak determined by Belsole et al. 
(200 4]). These a r e the s ame components denoted as Core and C 



in 



__ Burgett et a J ( 20041) . and as the NW and SW components in 
iDonnellv et alJ (120011) . Ou r components 1 and 2 correspond to 



a sparse component S by Hwang & Lee (2009), Belsole et al. 
(2004) shows an X-ray source in this direction. Hwang &Lee 
(2009) suggest that there has been past interacti on between the 
compo nents 1, 2, and 3. The component E in Hwang &Led 
(2009) corresponds approximately to our component 5, and may 
have interacted in the past with both the components 3 and 4. 
Hwang & Lee did not find significant differences between the 
galaxy colours in different components, similar to our study. 

Thus in the cluster A 1750 our analysis reveals a similar sub- 
structure to that found in other studies. 



4.3. The cluster A 1 773 in the outskirts of the supercluster 
SC1 126 

The next cluster in our sample is the X-ray cluster A1773, the 
member of the supercluster SCI 126, which is located in the out- 
skirts of the supercluster. The sky distribution of galaxies in this 
cluster is irregular, especially in the outer regions (Fig. [35] left 
panel). The 3D distribution of galaxies (Fig. [35] right panel) 
shows a component or tail of mostly elliptical galaxies and 
BRGs with negative peculiar velocities. The histograms of the 
peculiar velocities (Fig. [36b confirm that the brightest galaxies 
in this cluster have mostly negative peculiar velocities. They 
are located in one part of the sky distribution (see Fig. [35] left 
panel). The first ranked galaxy of the cluster is also located here, 
but has only a low peculiar velocity (v pec - -76 km/s). However, 
this subsystem of galaxies does not form a completely separate 
component; according Mclust, this cluster has one component 
with the mean uncertainty of the classification 0.041. The distri- 
bution of the peculiar velocities of spiral galaxies in this cluster 
is smoother than that of elliptical galaxies. 

The fraction of red galaxies in the A1773 is 0.67 only. The 
scatter of colours of red ellipticals is smaller than that of red spi- 
rals (with the rms values of 0.032 and 0.047, correspondingly, 
see Tableland Fig. [37). If galaxies with negative peculiar ve- 
locities have only recently joined this cluster they have probably 
finished their colour evolution before joining the main cluster. 

4.4. The cluster A1 809 in the outskirts of the supercluster 
SC1 126 

The last cluster in our sample is the cluster A 1809 in the out- 
skirts of the supercluster SCI 126. The sky distribution of galax- 
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Fig. 36. Distribution of the peculiar velocities of galaxies in the 
cluster A1773. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. The small dash in the upper panel indicates 
the peculiar velocity of the first ranked galaxy. 
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Fig. 37. Colour-magnitude diagram of galaxies in the cluster 
A 1773. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 
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Fig. 38. 2D and 3D views of the cluster A1809. Symbols are the 
same as in Fig. [3] 



ies in this cluster (Fig. [38] left panel) shows that this cluster con- 
sists of two components in the sky. Mclust also confirms that 
this cluster consists of two components with the mean uncer- 
tainty of classification 1.3 • 1CT 2 . In the main component of the 
cluster galaxies are concentrated in the center, and most BRGs 
of the cluster are also located here. The other component con- 
tains mostly loosely located spiral galaxies; there are also two 
elliptical galaxies and three BRGs, two of them located between 
this and the main component of the cluster. However, a word of 
caution is needed - in the poor component the number of galax- 
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Fig. 39. Distribution of the peculiar velocities of galaxies in the 
cluster A1809. Top down: all galaxies, elliptical galaxies, spiral 
galaxies, and BRGs. The small dash in the upper panel indicates 
the peculiar velocity of the first ranked galaxy. 
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Fig. 40. Colour-magnitude diagram of galaxies in the cluster 
A1809. Empty circles denote elliptical galaxies, crosses spiral 
galaxies, filled circles BRGs. 



ies is small, 17, the galaxies here are loosely distributed; it is 
possible that these galaxies cannot be considered as forming a 
reliable separate component. The 3D distribution of galaxies in 
the cluster A1809 (Fig. [38] right panel) and the distribution of 
the peculiar velocities of galaxies (Fig. [39]) show that in the main 
component of the cluster the galaxies have a bimodal distribution 
of the peculiar velocities, which is best seen in the distribution 
for spiral galaxies. The velocity histograms show that some el- 
liptical and spiral galaxies form a tail in the distribution of the 
peculiar velocities with high (v pec > 1000 km/s) values of the 
peculiar velocities. The distribution of the peculiar velocities of 
red spiral galaxies follows those of all spirals. The velocity dis- 
persion of the BRGs in this cluster is small, these galaxies are 
concentrated near the centre of the cluster. The distribution of 
the peculiar velocities of galaxies suggests that this cluster has 
formed by merging of two groups, and the third one (the poor 
component in the cluster A1809) is perhaps still infalling. 

The fraction of red galaxies in the cluster A1809 is 0.84, 
similar to that in the clusters in the core of the supercluster. 
The scatter of colours of red ellipticals and red spirals (with the 
rms values of 0.029 and 0.032, correspondingly, see Tableland 
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Fig. l40t is small. This suggests that the galaxies in the groups 
which formed the cluster A 1809 possessed their final colours 
before merging. 

The cluster A 1809 has been studied for subclustering by 
lOegerle & Hilll (11994 l200lh . who found no substructure here. 
The reason for the difference between their results and ours may 
be the difference in the data used - in their study the cluster in- 
cludes a smaller number of galaxies, and the smaller component 
is not well seen. As we mentioned above, it is also possible that 
the galaxies in the smaller component do not form a reliable sep- 
arate component. 



5. Discussion 

5.1. Dynamical state of rich clusters 

Our calculations with Mclust showed that among the ten clusters 
studied here five consist of more than one component, and the 
clusters A 1205 and A 1750 have the largest number of substruc- 
tures. Even one-component clusters have several components in 
the distribution of the peculiar velocities of galaxies. In one clus- 
ter, A1205, one of the components is mostly populated by red 
galaxies, both ellipticals and spirals. In other clusters, the distri- 
bution of the peculiar velocities of galaxies shows that elliptical 
galaxies are more concentrated towards the cluster centre or to- 
wards the centres of subclusters, than spiral galaxies. This is pos- 
sibly a signature of recent mergers that affects the morphology 
of galaxies in cluster. In the cluster A 1773 we see the opposite 
- the velocity dispersion of spiral galaxies is lower than that of 
elliptical galaxies. This shows that the formation history of our 
clusters was different. 

A word of caution is needed concerning the poorest clusters 
in our sample, the clusters A1650 and A1658: owing to the small 
number of galaxies in these clusters our results concerning them 
should be taken as a suggestion only. Even clusters with smaller 
numbers of galaxies have been studied for substruc tures (see, 
for examp le, Solan es et al.lfl9 99: Boschin et al. 2008); however, 
iBoschin et al.l (2008) mention that clusters with about 30 mem- 
ber galaxies are too small for the analysis of the substructure 
(see also t he discussion a bout the substructure statistics in small 
samples in Bivi ano et al.| [2006). 

Clusters with multiple components belong to all three su- 
perclusters of the SGW. Based on our earlier studies of the su- 
percluster SCI 126, we assumed that clusters in the high-density 
core of this supercluster are dynamically more evolved than clus- 
ters in other superclusters of the SGW, and close to be virialized. 
Now we see that this is not so, clusters in the high-density cores 
of superclusters also have substructure and thus are not yet viri- 
alized; 

Ragone-Fi gueroa & Plion is (2007) showed that in simula- 
tions group-size haloes located close to massive haloes have 
substructure. Plionis (2004) found that dynamically active clus- 
ters are more strongly clustered than the overall cluster popu- 
lation. Some studies of N-body models suggest that the frac- 
tion of haloes with substructure typically increas es in high- 
density regions (Ragone-Figueroa & Plionis! I2007L and refer- 
ences therein). In our study all clusters lie in superclusters, but 
we did not detect clear differences in the properties of clusters in 
the supercluster cores and in the outskirts or in the poor super- 
cluster SCI 91. 

We analysed the distribution of the peculiar velocities of 
galaxies in clusters. The number of galaxies in the clusters is 
not large, thus for the distribution of the peculiar velocities of 
galaxies we used the Shapiro- Wilk test to check the agreement 



with normality. In Table|2]we give the p-values of this test. These 
values show that the highest probability that the distribution of 
the peculiar velocities of galaxies in a cluster is not normal, is 
recorded for the cluster A 1066, one-component cluster, in which 
the distribution of the peculiar velocities of galaxies shows three 
components, a possible indication of merging subgroups. 

In Table|2]we give the values of kurtosis and skewness for the 
distributions of the peculi a r velocities of ga laxies i n clusters (see 
also [Solanes et al.l 119991: lOegerle & Hilll 1200 1[ IBoschin et al] 
120061: iBarrena et al.l 120071: iHwang & Led 120071) . The positive 
values of kurtosis in Table [2] suggest that all distributions are 
peaked. The skewness shows that the distributions are asymmet- 
rical, being skewed to the right or left, depending on the cluster. 
However, the /^-values are higher than 0.05 showing that devia- 
tions from the normal distribution are statistically not very sig- 
nificant. The /?-value for kurtosis for the cluster A 1205 is low, for 
skewness it is high. These tests compare different aspects of the 
distributions, thus the results may not be similar. We also see that 
the results of the tests alone are not always good indicators for 
the possible substructure in the distributions. We saw that in the 
clusters with several components in the sky the distributions of 
the peculiar velocities of galaxies in different components some- 
times (partly) coincide. Thus the low statistical significance for 
the deviations from normality is not surprising. 

Another indicator of virialization of clusters is the peculiar 
velocity of the first ranked galaxy in cluster. In six clusters from 
our sample the peculiar velocities of the first ranked galaxies 
v pec > 180 km/s. The relative peculiar velocities of the first 
ranked galaxies in our clusters, \v pec \/o- v (Table 13 can be di- 



vided into three classes: four clusters have \v pec \/o- v < 0.20, 
two - \v pec \l<j v - 0.41, and four clusters - \v pec \/cr v > 1.0. 
The low value of the relative peculiar velocity for the clusters 
A 1 205 and A 1 809 is misleading - the presence of multiple com- 
ponents in these clusters is a clear evidence of the nonvirialized 
state. Two other clusters with relative peculiar velocities less 
than 0.20<x,, are both one-component clusters. Both clusters with 
VpeJo~v — 0.41, although they are one-component systems, have 
several components in the distribution of the peculiar velocities 
of galaxies. Four multicomponent clusters have the highest val- 
ues of the relative peculiar velocities. 

High peculiar velocities of the first ranked g alaxies in clus- 
ters have been found in several recent studies (Oegerle & Hilll 
2001; ICoziol et al.l 120091) . Some of the cluster s in our sa mple 
and in the sample of rich clusters in ICoziol et al.1 (2009) co- 
incide; for common clusters our results agree well. The high 
peculiar velocities of the first ranked galaxies in clusters sug- 
gest that those clusters have not been yet virialized after a re - 
cent merger (JMalumuth 119921) . see also lOegerle & Hilll d2001l) ; 
ICoziol et a l. (2009). Before merger, these galaxies were the first 
ranked galaxies of one group which merged to form a larger sys- 
tem. As evidence for that, we found that in several of our groups 
the first ranked galaxy is located near the centre of one com- 
ponent. In this respect the cluster A 1205 is exceptional - in this 
cluster the first ranked galaxy lies between different components. 

All rich clusters we studied correspond to Abell clusters. An 
identification of clusters is usually complicated because of the 
differences between the data in the catalogues. This is also the 
case here. We compared the data of our clusters with the data of 
Abell clusters from the catalogue by Andernach & Tago (2009, 
private communication). This comparison showed that often sev- 
eral our groups can be associated with one Abell cluster. Our 
cluster, which we identified with the Abell cluster A 1066, is the 
richest cluster among the four groups which can be associated 
with this cluster; among the other groups one has 35 member 
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galaxies, others are pairs of galaxies. Other clusters from the cat- 
alogue by Andernach & Tago (2009) can be associated with one 
rich cluster and one to four poor groups with up to five member 
galaxies. This also confirms the presence of substructure in rich 
clusters. 

The cluster A 1750 has been studied for substructure by other 
authors, too. We showed above that the components found in 
this cluster by us coincide well with those found in other stud- 
ies. We found an especially good ag reement between our results 
and those by Hwang & Lee (2009), who used used a A - test 
(Dressier & Shectman 1988) to determine the substructure in the 
cluster A 1750. This shows that Mclust, probably applied here for 
the first time for the purpose to find substructure in galaxy clus- 
ters, gives results in agreement with other methods. 

The good agreement with other studies also supports our 
choice of the parameters for the friend-of-friend algorithm for 
group definition and suggests that in our catalogue groups with 
substructure are real groups and not complexes of small groups 
artificially linked together by an unreasonable choice of linking 
parameters. This example supports our opinion that in the T10 
the linking parameters were chosen reasonably. 

An increasing number of studies have recently 
shown the p resenc e of substructu r e in rich clusters 
(ISolanes et all 119991: lOegerle & Hilll 120011; iBurgett et al 



200 4 i Boschin et alJ 120061: 



Barrena et al. 2007; Hwang & Lee 



2007; Aguerri & Sanchez-Jansse nll2010h and in poor clusters 



Boschinetal. 2008); Vennik & Hoppf d2009l) found a sub 



Tovmassian & Plionis (2009) 



d 

structure also in poor groups. 

studied the properties of poor groups from the SDSS survey 

and showed that many groups of galaxies presently are not in 

a dynamical equilibrium but at various stages of virialization. 

Niemi et al. (2007) showed that a significant fraction of nearby 

groups of galaxies are not gravitationally-bound systems. This 

is important because the masses of observed groups are often 

estimated assuming that groups are bound. 

The presence of substructure in rich clusters, signs of possi- 
ble mergers, infall, or rotation suggests that the rich clusters in 
our sample are not yet virialized. The high frequency of these 
clusters tells us that mergers between groups and clusters are 
common - galaxy groups continue to grow. In high-density re- 
gions groups and clusters of galaxies form ear ly and could be 
more evolved dynamically (Temp el et al.ll2009l) . but here again 
the possibility of mergers is high, thus groups and clusters in 
high-density regions are still assembling. 



5.2. Galaxy populations in rich clusters 

The study of the galaxy content o f the rich clusters s howed in 
agreement with our earlier study (Ein asto et al.l 120081) that the 
fraction of red galaxies in clusters from the core of the super- 
cluster SCI 126 is very high, F rei ] > 0.8. However, this is as large 
as the fraction of red galaxies in the cluster A 1809 from the out- 
skirts of this supercluster, as well as in the cluster A1516 from 
the supercluster SCI 111. We plot the fractions of red galaxies 
in clusters and the fractions of elliptical galaxies in Fig. [41] This 
figure shows that the fraction of elliptical galaxies in clusters is 
proportional to the fraction of red galaxies. At the same time 
the luminosities of clusters vary strongly, i.e. the fraction of red 
galaxies in clusters from our sample does not correlate well with 
the luminosity of clusters. In our earlier studies we found that 
the fraction of red galaxies is the highest in clusters in the cores 
of rich superclusters , both for low- and high-luminosity clusters 
dEinasto et al.ll2007al) . 
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Fig. 41. Fractions of red galaxies and elliptical galaxies for our 
clusters. Symbol sizes are proportional to the luminosity of clus- 
ters. 



We found that approximately 1/3 of red gal axies are spi- 
rals. Recent results of the Galaxy ZOO project dMasters et al.l 
2009) showed that indeed a large fraction of red galaxies are 
spirals, and almost all massive galaxies are red independently of 
their morphological type. The colour-magnitude diagrams show 
that in some clusters (A1066, A1205, A1424) the scatter of the 
colours of red spirals is larger than the scatter of the colours of 
red ellipticals, suggesting that red ellipticals form a more homo- 
geneous population than red spirals. 

The number of BRGs in our clusters varies from four 
in A1658, which is the poorest cluster in our sample, to 
17 in the richest cluster, A1750. Some BRGs are spirals. 
BRGs can be found both in the central areas and in the out- 
skirts of clusters. The result that BRGs lie in clusters is con- 
sistent with t he strong small- scale clust e ring of the LRGs 
dZehavi et alJ 120051: lEisenstein et all 120051: iBlake et alJ [2008h 
that has been interpreted in the framework of the halo occu- 
pation model, where more massive haloes host a larger num- 
ber of LRGs dvan den Bosch et alJ 1 20081: IZheng et al.1 120081: 
iTinker et alJl2009UWatson et alj|201ob . 



The colours of galaxies in cluster environment may be af- 
fected by many factors. Colours may reach their observed val- 
ues during a merging event when galaxies are affected by pro- 
cesses in the cluster environment dBerrier et al.1 12009). For ex- 
ample, galaxies may be moved to the red sequence through the 
quenching of st ar formation in blue galaxies in the cluster envi- 
ronment ( Bow er et al.ll 19981: iRuhland et al]|2009t: iBamford et ail 
2009; Skelton et al. 2009). As evidence for that we found that in 
some of the clusters in our study the scatter of colours of galax- 
ies is large, probably due to merging of components that affects 
the colours of galaxies. 



At the same time ITinker et alJ (l2009h mention in the analy- 
sis of the distribution of distant red galaxies in the halo model 
framework that distant red galaxies must have formed their 
stars before they become satellites in haloes. This is possible if 
these galaxies were members of poorer groups before the groups 
merged into a richer one. We also find that in some clusters the 
range of red colours of galaxies is very small although these 
clusters consisted of several (probably merging) components. 
Probably in these cases the galaxies have finished their colour 
evolution in small groups before merging. 
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Table 4. The properties of the five most massive dark matter haloes 
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Columns in the Table are as follows: Col. 1: halo ID index, Col. 2: number of subhaloes (SH) at z = 0. Col. 3: mass of the main halo at z = (in 
units of 10 14 /r'M o ), Col. 4: the formation time, defined as the redshift when the halo has half of its mass at the present epoch, Col. 5: size of the 

halo (virial radius, in units of h' 1 Mpc), Col. 6: mean value and standard deviation of the formation time of subhaloes, Col. 7: total number of 
recent (z < 0.25) merger events ~LN merg , Col. 8: number of major mergers at z < 0.5, Col. 9: number of components in a halo at the present epoch 

found by Mclust as for observed clusters. 



5.3. Merger analysis of dark matter haloes from simulation 

5.3.1. Description of simulations 

To illustrate the formation of present-day haloes via multiple 
mergers during their evolution in a cosmological simulations we 
present in this subsection merger histories of dark matter haloes 
in an N-body simulation. 

We ran a ACDM simulation to study the merging history of 
dark matter haloes in a cosmological s i mulation. We u sed the 
GADGET-2 code (ISpringeletalJl200U ISpringelll2005h . In the 
simulation the dark matter haloes were i dentified with an al- 
gorithm called AHF (Amiga Halo Finder Knollmann & Knebe 
2009), which is based on the adaptive grid structure of the simu- 
lation code in AMIGA (Adaptive Mesh Investigations of Galaxy 
Assembly). All haloes consist of the main halo - the most mas- 
sive halo in the group, and subhaloes located within the virial 
radius sphere of the main halo. 

For the simulation we adopted a cosmological ACDM model 
with (Q. m + £2a + Q.I, — 1) and h - 0.71, the dark matter density 
Qdm - 0.198, the baryonic density Q./, - 0.042, the vacuum 
energy density Q.^ - 0.76 and the rms mass density fluctua- 
tion parameter <x 8 = 0.77. The simulation has 512 3 dark mat- 
ter particles, the volume of the simulation is 80/z~' Mpc 3 . In 
this simulation the masses of the largest haloes are of the or- 
der of 10 l *h~ l M Q , the mass of an individual DM particle is 
2.54xlO 8 /z~ 1 M . Our requirement for halo identification is that 
it has at least 100 particles, and therefore the minimum mass 
for haloes is 2.5xl0 I0 /i~ l M Q . In the simulation, the main haloes 
with embedded subhaloes correspond to the observed groups and 
clusters of galaxies. Studies of substructure in observed groups 
and clusters of galaxies and in dark matter haloes in a similar 
manner provide an important probe for the formation of ob- 
served galaxy systems. 



5.3.2. Halo mergers 

The hierachical formation theory predicts that small haloes form 
first and large haloes are later assembled via mergers. To study 
the mergers between haloes during their formation, we calcu- 
lated complete merger histories for the five richest (most mas- 
sive) main haloes in the simulation. We focused on the late-time 
evolution, because it is more likely that the events occuring dur- 
ing this epoch may still be seen in the structure of the haloes. 

In Table |4] we show the results of our analysis. We calcu- 
lated the formation times for the main haloes, defined as the red- 
shift when halo has half of its mass at the present epoch, and the 
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Fig. 42. Complete merger tree for a halo with ID = 1. Major 
merger events are indicated by circles with the size proportional 
to the mass of the merged halo. 



mean value of the formation time for all the subhaloes with the 
standard deviation. We also show in Table|4]the total number of 
merger events for main haloes for z < 0.25 and the number of 
major mergers where the mass of the merging halo is at least 0.1 
of the main halo mass. 

To illustrate the complexity of halo histories, we show two 
examples of merger trees that we used for calculating the prop- 
erties of mergers. In Fig. |42]we show the merger tree for the halo 
ID 1, which had three large recent merger events (at redshifts z = 
0.0, 0.049, and 0.211). This halo has recently accreted its mass 
and its formation time is z — 0.271. This can be compared with 
the merger history for the halo ID 3, which has the same forma- 
tion time and almost the same amount of recent mergers, but in 
this halo only one large merger event at redshift z — 0. 1 is clearly 
visible. Thus the merger tree of halo ID 3 is quite different from 
that of halo ID 1 (Fig. l43l . Interestingly, the mean value of the 
formation time for subhaloes in the halo ID 3 is lower than that 
of the halo ID 1 . The number of components found by Mclust is 
2 for the halo ID 1 and 1 for the halo ID 3 (Tableland Fig. Eg). 
Qualitatively one can say that more active late-time merging is 
visible in the clustering properties of haloes at z — 0. 
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Fig. 43. Complete merger tree for a halo with ID = 3. Merger 
events as indicated in Fig. [42] 
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Fig. 44. Mclust results for haloes 1 and 3: 3D view of the com- 
ponents. 

A similar analysis was carried out for the five most massive 
simulated haloes. Our results show that a late-time formation 
of the main haloes and a number of recent major mergers can 
cause a late-time subgrouping of haloes. A large scatter in the 
formation times for subhaloes hints towards late time merging 
(the standard deviation in the Col. 6). In general, we notice that 
the number of recent large merger events is the best indicator for 
additional clustering and that it is strongly correlated with the 
Mclust findings. To better understand the merging histories of 
dark matter haloes a study of a larger sample of haloes is needed. 

6. Conclusions 

We studied the properties of rich clusters in different superclus- 
ters in the SGW. Our conclusions are as follows. 

1) We showed with Mclust, a publicly available package in the 
R statistical environment, that the richest clusters in all su- 
perclusters of the SGW have substructure. 

2) There are components in the distribution of the peculiar ve- 
locities of galaxies in clusters both with multiple components 
and in one-component clusters - evidence of merging com- 
ponents. 

3) The peculiar velocities of the first ranked galaxies in clusters 
with multiple components are high. 



4) Most rich clusters in the supercluster SCI 126 and the richest 
cluster in the supercluster SCI 111 have a very high fraction 
of red galaxies (F r > 80 %). About 1/3 of red galaxies in 
clusters are spirals. In several clusters the scatter of colours 
of red elliptical galaxies in the colour-magnitude diagram is 
smaller than the scatter of colours of red spirals, suggesting 
that red elliptical galaxies form a more homogeneous popu- 
lation than red spiral galaxies. The presence of colour-related 
substructure in these clusters suggests that the colour evolu- 
tion of galaxies in these clusters has already finished before 
merging into clusters. 

5) Our clusters host a large number of BRGs. Some BRGs are 
spirals. They can be found both in the central areas and in 
the outskirts of clusters. 

6) We calculated the complete merger trees for the richest dark 
matter haloes in an N-body simulation and showed that the 
late-time formation of the main haloes and the amount of 
recent major mergers can cause the late-time subgrouping of 
haloes. 

Our results about the presence of substructure in rich clus- 
ters, signs of possible mergers, infall or rotation suggest that 
rich clusters in different superclusters of the Sloan Great Wall, 
including the cores of superclusters, are not yet virialized. The 
richest clusters in the SGW are still assembling. 

As a next step in the study of substructures in groups and 
clusters, we plan to study a larger sample of groups and clus- 
ters, including poor groups of galaxies and low-density global 
environments. 
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